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We propose the local density approximation (LDA)+Gutzwiller method incorporating Green's 
function scheme to study the topological physics of correlated materials from the first-principles. 
Applying this method to typical mixed valence materials SmBg , we found its non-trivial Z2 topology, 
indicating that SmBe is a strongly correlated topological insulator (TI). The unique feature of this 
compound is that its surface states contain three Dirac cones in contrast to most known TIs. 

PACS numbers: 



Most of Z2 topological insulators (TI) discovered 
up to now are semiconductors which are free of strong 
correlation effects and their topological nature can thus 
be predicted quite reliably by first principle calculations 
based on density functional theory (DFT)@. The Z2 
classification of band insulators has been generalized to 
interacting systems by looking at its response to external 
electric magnetic field, namely the topological magneto- 
electric effect (TME)0|. A correlated insulator is a TI 
if the 6>-angle defines in TME is n. Given the TME the- 
ory and several simplifications however, its applica- 
tion to realistic materials is still absent due to: (1) the 
lack of suitable compounds; (2) the difficulty to compute 
reliably correlated electronic structures from the first- 
principles. In this letter, we study a special class of mate- 
rials, the mixed valence (MV) compounds, which contain 
rare-earth elements with non-integer chemical valence. 
By combining the Gutzwillcr variational approach from 
the first-principles and the Green's function method for 
the TME, we found SmBg is a 3D correlated TI. Interest- 
ingly, it has three Dirac cones on the surface, in contrast 
to most of the known TIs. 

Classic MV compounds, such as SmBg and YbB^Q, 
share the following common features, (i) The x-ray 
photoelectron spectroscopy (XPS) and X-ray absorption 
spectra (XAS) contain peaks from both divalent and 
trivalent multiplets with comparable spectral weig ht[l3- 
[l2|, indicating the valence of Sm or Yb to be close to 2.5. 
(ii) A small semiconducting gap opens at least in part of 
the Brillouin zone (BZ) at low temperature 13]. Previous 
electronic structure studies indicate that the electrons 
transfer from Sm (or Yb) 4/ orbitals to 5d orbitals, which 
leads to fractional occupation in 4/ orbitals [3, EH ■ From 
the band structure point of view, the electron-transfer be- 
tween 4/ and 5al orbitals indicates possible band inversion 
between them, which is a crucial ingredient to realize TI 
(if it happens odd number of times in the entire BZ). 
At the mean while, SmB 6 has been recently suggested 



(A) 







Xi 


M 


/ 








R 








r 








K 


V 

A 




M 













FIG. 1: (A) The CsCl-type structure of SmB 6 with Pm3m 
space group. Sm ions and B6 octahedron are located at the 
Conner and center of the cubic lattice respectively. (B) The 
bulk and surface BZ for SmB6. 



to be among the possible candidates of the interesting 
topological Kondo insulators flol - flll ] . Although the ana- 
lytical studies of MV compounds have been conducted by 
several groups, the reliable first principle studies, which 
is crucial to identify the possible topological phases, is 
still lacking, due to the strong correlation nature of these 
materials flfA [20|. Here we report that the local density 
approximation (LDA)+Gutzwiller method, a newly de- 
veloped first-principles tool for correlated electron sys- 
tems, enables us to search for the topological phases in 
correlated matters. By taking the MV compound SmB 6 
as an example, we will focus on two key issues: (i) how 
to compute Z 2 topological index with LDA+Gutzwiller; 
(ii) what's the topological nature of SmBg with the strong 
correlation effects among the / electrons? 

The LDA+Gutzwiller method combines the DFT 
within LDA with the Gutzwiller type trial wave func- 
tion, which takes care of the strong atomic features in 
the ground state. Here wc just sketch the most impor- 
tant aspects, and leaving details to reference [2l"1423l | . We 
start from the common Hamiltonian (used for most of 
the LDA++ schemes), 
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where H^oa is the single particle Hamiltonian obtained 
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by LDA and Hint is the local interaction term for the 
4/ electrons described by the Slater integrals Fq, F 2 , F±, 
F 6 . In the present study, for Sm atom we choose Fo=5.8 
eV and F 2 = 9.89eV, F = 7.08eV,F 6 = 4.99eV to be 
their atomic values |2J, |25[. Hoc is the double count- 
ing term representing the interaction energy already con- 
sidered at the LDA level. In the present paper, we 
compute the double counting energy using the scheme 
described in reference |26j|. We then use the following 
Gutzwiller trial wave function, \G) = P G \0) = ^ ^1°). 
where |0) is a non-interacting state (obtained from LDA), 
Pi = J^ r . Arjr^I^I is the Gutzwiller projector at the 
i-th site with being the atomic eigenstates and Ar be- 
ing the variational parameters to be determined by min- 
imizing the ground state total energy (under Gutzwiller 
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LDA+Gutzwiller 



approximation) [22t l27l |28| , E G =(0\H eff \0) + J2 r x rE r . 



Here H e ff = PgHldaPg is called renormalized effec- 
tive single particle Hamiltonian for the quasi-particles 
and Ey is the eigen energy of the T-th atomic eigenstate. 
The scheme preserves the nice aspect of being variational, 
however, it is beyond LDA (and also LDA+C/) because 
the total energy Eq now relies on the balance between 
the renormalized kinetic energy of quasi-particlc motion 
and the local interaction energy, which is configuration- 
dependent. We note that including the complete form of 
all the Slater integrals (F to Fq) in the local interaction 
is essential to obtain the correct electronic structure for 
these materials. 

The linear response theory for the coefficient of TME 
has been developed and simplified by Z. Wang et al^. 
For interacting systems, when the self-energy contains no 
singularity along the imaginary axis, the formula for the 
TME coefficient only requires the single particle Green's 
function at zero frequency, g{k, 0). Because the singu- 
lar point for the self-energy along the imaginary axis 
only appears for a Mott insulator with completely lo- 
calized /-orbitals, which is not the case for SmBg, we 
can safely apply the above method here. Therefore the 
way to determine the Z 2 invariance is just diagonalizing 
the Hermite matrix — g(/c,0) _1 = Ho — /i/ + S(0) and 
treating the eigenstates with the negative eigenvalues as 
the "occupied states" . If the system also has spacial in- 
version center, the Z 2 invariance can be simply deter- 
mined by counting the parities of these "occupied state" 
on the time reversal invariant momenta (TRIM) points, 
the same as that done for the non-interacting topological 
band insulators [H, 0|. 

In the Gutzwiller approximation (or the equivalent 
slave boson mean field approach), the low energy single 
particle Green's function of an interacting system can be 
expressed by the quasi-particle effective Hamiltonian and 
the quasi-particle weight z (which in general is a matrix) 
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FIG. 2: The calculated band structure for SmBe in different 
energy scales by LDA (A, B) and LDA+Gutzwiller(C, D).The 
B and D are just the zoom in of A and B respectively, around 
the fermi level. Compared to the LDA results, the 4/ j=7/2 
bands in LDA+Gutzwiller have been pushed up to be about 
4.0 eV above the fermi level, leaving the band structure near 
the fermi level being dominated by 4/ j=5/2 and 5d states. 
The blue, red and grey colors represent the weight of the 
orbital character for 5d, 4/ and 2p respectively. 



where the second term describes the incoherent part 
of the Green's function, which is ignorable for low fre- 
quency. Therefore by applying the green's function 
method described above, we reach the conclusion that 
the Z 2 invariance in LDA+Gutzwiller method is deter- 
mined by the occupied eigenstates of the Gutzwiller ef- 
fective Hamiltonian H e ff, which can be interpreted as 
the band structure of the "quasi-particles" . 

We now focus on SmBg, which is crystallized in the 
CsCl-type structure with Sm ions and Bq clusters be- 
ing located at the Conner and body center of the cubic 
lattice respectively (see Fig.l). The LDA part of the 
calculations has been done by full potential linearized 
augmented plane wave (LAPW) method implemented in 
the WIEN2k package [29j. BZ integration was performed 
on a regular mesh of 12x12x12 k points. The muffin-tin 
radii(i?MT) of 2.50 and 1.65 bohr were chosen for Sm 
and B atoms, respectively. The largest plane-wave vec- 
tor K max was give by RuT^-max = 8.5. The spin-orbit 
coupling (SOC) is included in all calculations. 

The band structure obtained by LDA is shown in FigJ2] 
(A and B), where we can find three major features. 

(i) The Sm-4/ orbitals, which split into the j=5/2 and 
j=7 1 1 manifolds due to the SOC, form narrow bands re- 
spectively with width around 0.5 eV near the Fermi level. 

(ii) The low energy band structure is semiconductor-like 
with a minimum gap about 15meV along the T-X direc- 
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tion. (iii) There is clear band inversion features at the 
X-points, where one 5d band goes below the / bands, 
which reduces the occupation number of the /-states to 
be around 5.5. 

The first two features of the LDA band structures 
are not consistent with the experimental observations. 
Firstly, the XPS measurements find quite strong 4/ 
multiplet peaks, indicating strong atomic nature of 4/ 
electrons in SmB 6 (in other words, most of the 4/ 
electrons are not envolved in the formation of energy 
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bands) [1 
optical [3 




and 



Secondly, both transport 
361 ] measurements reveal the formation of a 
small gap only for temperature below 50 K. The above 
two features imply that the / electrons in SmBg have 
both localized and itinerant natures and the correct de- 
scription of its electronic structure should include both 
of them. The qualitative physical picture can be ascribed 
to Kondo physics, which involves both itinerant and lo- 
calized electronic states [13] ■ At high temperature, these 
localized orbitals (4/ states) are completely decoupled 
from the itinerant energy bands, forming atomic multi- 
plet states. While in the low temperature, the coherent 
hybridization between localized and itinerant states is 
gradually developed leading to the formation of "heavy 
quasi-particle" bands. The insulating behavior appears 
when the chemical potential falls into the "hybridization 
gap" between the heavy quasi-particle and the conduc- 
tion bands, which is mainly of the 5d orbital character in 
SmBg. 

Such Kondo picture can be nicely captured by 
our LDA+Gutzwillcr calculations, which provide equal- 
footing descriptions of both the itinerant 5d bands 
and those heavy quasi-particle states formed by 4/ or- 
bitals. The quasi-particle band structures obtained by 
LDA+Gutzwillcr is shown in Figj2](C and D). Compared 
with the LDA results, there are three major differences 
induced by the strong correlation effects, (i) The 4/ 
j=7/2 bands are pushed up to be around 4.0 eV above 
the fermi level, leaving the band structure near the fermi 
level being dominated by 4/ j=5/2 and 5d bands, (ii) 
The band width of those "heavy quasi-particle" bands 
formed by 4/ orbitals are much reduced to be less than 
0.1 eV. The quasi-particle weight z of these "heavy quasi- 
particle" bands is about 0.18, indicating that the quasi- 
particles are formed by less than 20% of the 4/ spec- 
tral weight and the remaining weight is attributed to the 
atomic multiplets (or the Hubbard bands), (iii) A hy- 
bridization gap between 4/ quasi-particle bands and the 
itinerant 5c? bands appears along the F to X direction. 
The detail hybridization process is illustrated in Fig. [3] 
Along the r to X direction, the point group symmetry 
is lowered to be C^ v , whose double group contains two 
irreducible two-dimensional representations, Tg and Ty. 
Along r to X, the original 4/ j=5/2 orbitals split into 
two T7 and one Tq bands, which cross with another T-? 
band formed by 5d orbitals. The hybridization terms 
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FIG. 3: Band hybridizations between 4/ and 5d along V to 
X direction. The original 4/ j=5/2 orbitals split into two r£ 
and one Tg bands. The Tj bands cross and hybridize with 
another V% band formed by 5d orbitals. 



among the three Tr bands are allowed and generates the 
hybridization gap, which is around 10 mcV now. We note 
that in contrast to the LDA results, the semiconductor 
gap obtained by LDA+Gutzwiller is indirect, which is 
quite consistent with the transport measurements [36| . 

We would emphasize that by LDA+Gutzwiller, we can 
only obtain the quasi-particle part in the spectral func- 
tion, but not the Hubbard bands containing the atomic 
multiplet structure [H, [2(|. While the Gutzwiller type 
wave function can well capture the multiplet features 
in the ground state 22 , 27 , |3tJ • The Gutzwiller varia- 
tional parameter Ar determines nothing but the proba- 
bility of each atomic configuration T in the ground state 
(which is defined as (G|r)(r|G)). In Fig. H we plot 
the corresponding probabilities for SmBg obtained by 
LDA+Gutzwillcr together with that obtained by LDA 
wave function((0|r)(r|0)). The LDA ground state is 
dominated by the atomic configurations with the number 
of /-electrons Nf=6, on the other hand however, the dis- 
tribution of the probability of atomic states obtained by 
LDA+Gutzwiller is almost equally concentrated on two 
atomic multiplet states with five and six /-electrons re- 
spectively, leading to approximately +2.5 valence of Sm. 

Moreover, the Gutzwiller wave function provides cor- 
rect description of the intermixing between j = 7/2 and 
j = 5/2 orbitals, which can not be captured by LDA 
only calculation and manifest itself in the average occu- 
pancy of the f orbitals. The occupancy of j = 5/2 and 
j = 7/2 orbitals arc 5.31 and 0.22 respectively using the 
LDA type wave function (|0)). While the f orbital oc- 
cupancy is modified to be 3.64 for j = 5/2 and 1.89 for 
j = 7/2 orbitals using the Gutzwiller type wave func- 
tion (|G)). The dramatic increment of the occupancy 
for j = 7/2 orbitals is the important consequence of the 
Fi — Fq terms in the atomic interactions, which can not 
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FIG. 4: The probability of atomic eigenstates in the ground 
state obtained by (a) LDA alone and (b) LDA+Gutzwiller. 
Nf is the total number of / electrons for the corresponding 
atomic eigenstates. The horizontal axis denotes the corre- 
sponding atomic eigen energy. 

be expressed by a pure "density-density" form in any sin- 
gle particle basis and generate strong multi-configuration 
nature for the ground state wave function. 

TABLE I: The products of parity eigenvalues of the occupied 
states for TRIM points, E, X, R and M in the BZ. 

E 3X II 3M_ 

+ + + 



tight binding Hamiltonian and the same rotational in- 
variant Gutzwillcr approach on a 40-laycrs slab. The 
obtained quasi-particlc bands of the slab arc plotted in 
FiglU It is clearly seen that the surface states contain 
three Dirac cones: one cone is located at f and the other 
two are at two X points of the surface BZ. The inter- 
esting multi-Dirac-cones behavior is quite unique among 
the existing TIs and is a natural consequence of the band 
inversion at the bulk X points, which are projected onto 
to f and two X points of the (001) surface BZ. The 
multiple Dirac cones on the surface of SmBg may gener- 
ate interesting physical phenomcnas, such as the unique 
quasi-particle interference pattern in scanning tunneling 
microscope (STM), which will be studied in our further 
publications. 
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FIG. 5: The surface states of SmB6 on the (001) surface. 
The surface states are obtained by LDA+Gutzwiller with a 
40 layer slab on the basis of projected Wannier functions. 
(Inset) The fermi surfaces for SS on (001) surface. 



The band inversion feature around the X points is well 
persisted in the LDA+Gutzwiller quasi-particle bands. 
As discussed in previous paragraphs, the topological na- 
ture of this interacting system is fully determined by the 
Gutzwillcr effective Hamiltonian H e ff. Since the spacial 
inversion symmetry is present for SmBg , we are now able 
to determine its topological nature by simply counting 
the parities of those occupied quasi-particle states at 8 
TRIM points. As listed in TableU the parities are all pos- 
itive except the X points. Because there are totally three 
equivalent X points in the whole BZ, the topological 
index for SmBg has to be odd, resulting in a strongly 
correlated topological insulators with topological indices 
(Till). 

To see the topological surface states, we construct 
a tight bin ding model using the projected Wannier 
functions [38|-|40(, which can reproduce the LDA band dis- 
persion quite precisely. The surface states (SS) of the 
(001) surface are then obtained by combining the above 



In summary, we have developed the LDA+Gutzwiller 
method incorporating the Green's function scheme to 
study the the topological phases of strongly correlated 
materials from the first-principles (beyond LDA and 
LDA+U). This method is systematically applicable to 
all correlated compounds as long as the quasi-particlc 
weight is not reaching zero. Both quasi-particlc bands 
and atomic multiplcte structures can be well captured in 
the present technique. Applying this method onto typical 
mixed valence compound SmBg, wc demonstrate that it 
is a strongly correlated 3D TI with unique surface states 
containing three Dirac cones on the (001) surface. The 
strong interaction among the / electrons reduces both 
the band width and the quasiparticle weight for almost 
one order, but the topological feature remains. 
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